NN Scattering Phase Shifts in a Chiral Constituent Quark Model 



D. Bartz* and Fl. Stancu + 
Universite de Liege, Institut de Physique B.5, Sari Tilman, B-4000 Liege 1, Belgium 

(February 1, 2008) 

O" 

§ : Abstract 
(N 
O 

QJ ■ We study the nucleon-nucleon interaction within a chiral constituent quark 

Q 

£Sj ■ and 1 S'o phase shifts by using the resonating group method. They clearly 

>. 

indicate the presence of a strong repulsive interaction at short distance, due 



model which reproduces succesfully the baryon spectra. We calculate the 3 S" 



to the spin-flavour symmetry of the quark-quark interaction and of the quark 



o 

O 

■ interchange between the two interacting nucleons. A cr-meson exchange quark- 

o . 

quark interaction, providing a medium-range attraction, helps to get closer 
to the experimental phase shifts. 

CD ' 

PACS number(s): 24.85. +p, 21.30.-x, 13.75. Cs 

"X' 

Keywords : chiral constituent quark model, pseudoscalar meson exchange interaction, 

d ; 

nucleon-nucleon phase shifts, resonating group method 
Introduction 

The aim of the present work is to study the nucleon-nucleon (NN) interaction and in 
particular the NN scattering phase shifts in the framework of a chiral constituent quark 
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model which has proved succesfull in baryon spectroscopy In the calculation of the 

NN phase shifts we employ the resonating group method (RGM) |||7|] which has been first 
applied to nuclear physics. This method is very convenient for treating the interaction 
between composite particles. In the NN problem each nucleon is supposed to be composed 
of three quarks and the interaction between quarks is here provided by the same chiral quark 
model, both for quarks belonging to the same nucleon or to different nucleons. 

The study of the NN interaction in the framework of quark models has already some 
history. We restrict our discussion to nonrelativistic models where the RGM can be applied. 
Twenty years ago Oka and Yazaki || published the first nucleon-nucleon L = phase shifts 
calculated with the resonating group method. The work || and subsequent developments 
are revised in ||. Those results were obtained from models based on one-gluon exchange 
(OGE) interaction between quarks. Based on such models one could explain the short-range 
repulsion of the NN interaction potential as due to the chromomagnetic spin-spin interaction, 
combined with quark interchanges between 3q clusters. In order to describe the data, long- 
and medium-range interactions were added at the nucleon level. 

Here we employ a constituent quark model where the short-range quark-quark interaction 
is entirely due to pseudoscalar meson exchange instead of one-gluon exchange. This is the 
chiral constituent quark model proposed in Ref. Q and parametrized in a nonrelativistic 
version in Refs. [@||. A semirelativistic version is also presented in Ref. [|J. The present 
status of this model is summarized in Ref. ||. The spin-flavour symmetry structure of the 
model is getting support from the phenomenological analysis of L = 1 negative parity baryon 
TOf] . Also lattice |1]J and 1/Nc QCD studies Jl2| have a consistent interpretation 



resonances 



in a constituent quark model with pseudoscalar meson exchange interaction. 

The origin of the model is thought to lie in the spontaneous breaking of chiral 

symmetry in QCD which implies the existence of constituent quarks with a dynamical mass 
and Goldstone bosons (pseudoscalar mesons). According to the two-scale picture of Manohar 
and Georgi [T3L at a distance beyond that of spontaneous chiral symmetry breaking, but 



within that of the confinement scale, the appropriate degrees of freedom should be the 



constituent quarks and the chiral meson fields. If a quark-pseudoscalar meson coupling is 
assumed, in a nonrelativistic limit one obtains a quark- meson vertex proportional to a ■ q \ F 
with a the Pauli matrices, q the momentum of the meson and \ F the Gell-Mann flavour 
matrices. This generates a pseudoscalar meson exchange interaction between quarks which 
is spin and flavour dependent. In the following, the model ]!]-§[, based entirely on this 
interaction, is referred to as the Goldstone boson exchange (GBE) model. 

It is important to correctly describe both the baryon spectra and the baryon-baryon 
interaction with the same model. The model |]-[|] gives a good description of the baryon 
spectra and in particular the correct order of positive and negative parity states, both in 
nonstrange and strange baryons, in contrast to the OGE model. In fact the pseudoscalar 
exchange interaction has two parts : a repulsive Yukawa potential tail and an attractive 
contact ^-interaction. When regularized , the latter generates the short-range part of the 
quark-quark interaction. This dominates over the Yukawa part in the description of baryon 
spectra. The whole interaction contains the main ingredients required in the calculation of 
the NN potential, and it is thus natural to study the NN problem within the GBE model. 
In addition, the two-meson exchange interaction between constituent quarks reinforces the 
effect of the flavour-spin part of the one-meson exchange and also provides a contribution 
of a cr-meson exchange type required to describe the middle-range attraction. 

Preliminary studies with this interaction have been made in Refs. [p~5^p~7f| . This work is 



a natural extension of these previous studies. Ref. |T5| was rather exploratory about the 
role of a spin-flavour dependent interaction in giving rise to a repulsive core. Within the 
parametrization of the GBE model it was found that at zero-separation between two 3q 
clusters the height of the repulsive core is 0.830 GeV and 1.356 GeV in the 3 Si and 1 So 
channels respectively. The spin-flavour symmetry and the parametrization M of the GBE 



model favours the | [42]o[51]^,s) state, which becomes highly dominant. In Ref. |I6" |, instead 
of cluster model states, a better basis obtained from single-particle molecular type states, 
has been employed. The situation has been found to be similar, the repulsion being reduced 
by about 200 MeV in the 3 Si channel and by about 400 MeV in the l So channel. This is 



natural because the molecular orbital basis gives a lower bound of the expectation value of 



the Hamiltonian in the six-quark basis. In Ref. |I7j an adiabatic NN potential was calculated 
based on the model 0. It was found that none of the bases, cluster or molecular, leads to 
an attractive pocket. An attraction was simulated by introducing a a- meson exchange of an 
analytic form similar to that of the pseudoscalar meson exchange of [Q. 

To our knowledge, the present study is the first derivation of NN scattering phase shifts 
with the full SU(3) version of the GBE model which provides a comparison with the ex- 
perimental data. Similar RGM calculations 0], where a simplified SU(2) version of the 
GBE model has been used, have also shown the presence of a short-range repulsion in the 
behaviour of the phase shifts in the L = channel. However that paper did not aim at a 
comparison with the experiment. 

The model 

The GBE Hamiltonian considered below has the form |J 

9 

ff = E"N + £|!r-*G + £ Vconffa) + £ V x ( rij ) , (1) 
i i=l 1 i<j i<j 

where Kq is the kinetic energy of the center of mass. The linear confining interaction is 

Vcanffa) = -|a? ■ A • {Cr i3 + V ) (2) 
and the spin-spin component of the GBE interaction in its SUp(3) form is 



Vxira) = {£ V^XfXf + £ V K ( nj )\?\f + V^r^X* + ^(r y )}di • d$ . (3) 

F=l F=4 6 



The interaction (3) contains 7 = 7T,K,r), and rj' meson-exchange terms and Vyfcj) is given 
as the sum of two distinct contributions: a Yukawa-type potential containing the mass of 
the exchanged meson and a short-range contribution of opposite sign, the role of which is 
crucial in baryon spectroscopy. For a given meson 7, the exchange potential is 

a 2 1 p-th T p-^r 

= — a;^}, (4) 

4:71 lzminij 1 r 1 r 
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where A 7 = Ao + k/i 7 . For a system of u and d quarks only, as is the case here, the K 
exchange does not contribute. In the calculations below we use the parameters of Refs. ||. 
These are 

m u4 = 340 MeV, C = 0.77 far 2 , 
fi* = 139 MeV, fir, = 547 MeV, /v = 958 MeV, 

2 2 2 

9 JI1 = 9 M = L24 9jU = 2J652 

47T 47T 47T 

A = 5.82 fm -1 , k = 1.34, VJ, = -112 MeV . (5) 



The reason of using the parametrization |3|, instead of 0, as in the previous work fl5H17|, 
is that it is more realistic. Its volume integral, i.e. its Fourier transform at q = 0, vanishes, 
consistently with the quark-pseudoscalar meson vertex proportional to a • q \ F . In addition 
this interaction does not enhance the quark-quark matrix elements containing lp relative 
motion, as it is the case with the parametrization [f2[. This point has been raised in Ref. 
p~9|] . As a net result, in the parametrization [fj| the attraction due to the Yukawa-potential 
tail is overwhelmed by the repulsion resulting from the short-range part of the hyperfine 
interaction fl7 |. 



Some more comments are in order before proceeding further. 1) The above parametriza- 
tion gives a good description of baryon spectra. We do not change any parameter obtained 
from the fit 0. Such a parametrization is, of course, only effective. However, irrespec- 
tive of the parametrization, the flavour-spin symmetry is essential in this model. 2) There 
are also semirelativistic versions of the GBE model, as for example @] but first one has to 
study whether or not the resonating group method is applicable to semirelativistic six-quark 
Hamiltonians. The present RGM approach lies heavily on an s 3 structure of the nucleon 
wave function and such a simple description is inadequate for semirelativistic wave functions. 
Anyhow, for a fair comparison with results based so far on nonrelativistic OGE models a 
nonrelativistic version of the GBE model is desirable. 
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The nucleon 

In the RGM approach the wave function of the ground state nucleon must be known (see 
next section). As indicated above, its orbital part </> has an s 3 structure. This is a very good 
approximation to the exact wave function. The function <ft is fully symmetric with respect 
to any permutation of S3 and is chosen of the form 

</> = fl9(ri,b) , (6) 

i=l 

with g{fi, b) given by 

9(r,b) = {^f /A e~i , (7) 
nb z 

The size parameter b appearing in (0) is obtained variationally from the stability condition 
(see for example Ref. 0) 

^(WM = o, (8) 

where H is the Hamiltonian ([[J) written for a 3q system. This condition gives b = 0.44 fm, 
which we use below. 

The six-quark state 

From symmetry considerations (see for example ||20|| ) we can find which channels con- 
tribute to a totally antisymmetric six-quark state of a given spin S and isospin / when the 
orbital, spin, flavour and colour degrees of freedom of the nucleon are taken into account. 
For example for SI = (10) or (01) the channels are NN, AA and the hidden colour CC. In 
this work the discussion is restricted to NN channels only. Preliminary studies with coupled 
channels indicate that the influence of AA and CC channels on the phase shifts is very small 

PH. 
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The resonating group method 

Here we shortly describe the resonating group method. A detailed account of its appli- 
cation to the NN scatering can be found in the original papers by Oka and Yazaki ||. The 
resonating group method is a well established method for studying the interaction between 
two composite particles. First it has been succesfully applied to the derivation nucleus- 
nucleus interaction (see e.g. 0). Since the work of Oka and Yazaki, it is being used to 
calculate NN phase shifts assuming a quark structure for the nucleons. In the NN system 
each nucleon is treated as a 3q cluster. In the one channel approximation of the RGM the 
NN system is described by a 6q wave function of the form 

${U, g B , Rab) = A[<t>(g A )<t>(g B )x(RAB)} , (9) 

where £a = (^1,^2) and £b = (£3, £4) are the internal coordinates of the clusters A & B and 
Rab = Ra — Rb is the relative coordinate between the two clusters. Thus 4>i (i — A B) 
are the internal wave functions of the clusters which are supposed to be known and \ is the 
unknown wave function describing their relative motion. The development below is based 
on the assumption that <fti has a simple s 3 structure which is reasonable for nonrelativistic 
models. The antisymmetrization operator A is 

^=l-EE^r (10) 

t=l j=4 

From the variational principle one can obtain the equation determining the relative wave 
function x(Rab) as 

/ ^ + (U)4> + ^b)(H - E)^ A ^B,RAB)d 3 Ud^B = , (11) 

where H is the Hamiltonian ([]]) of the six-quark system. We introduce the Hamiltonian 
kernel 

H(R,R') = J <j) + (U)<t> + (tB)S(R- ^HA^iUm^SiR' ~ R A B)}d 3 Ud^Bd 3 R A B (12) 
and the normalization kernel 
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N{R,R>) = J <P + {U)<P + {£ B )KR- RAB)A[<P(i A )(f>(iB)S(R'-RAB)]d 3 Ud^Bd 3 RAB ■ (13) 
Then Eq. ( pi]) can be written as 

J C(R,R') X (R')d 3 R' = , (14) 

where C{R,R!) = TC(R, R!) - EAf(R, R'). This is the RGM equation. Here it has been 
solved by using the method of Ref. [[I]]. Accordingly, the relative wave function x{R) nas 
been expanded over a finite number of equally displaced Gaussians. For scattering states 
this expansion holds up to a finite distance R = R c , where R c depends on the range of 
the interaction. Beyond R c , x{R) is written as the usual combination of Hankel functions 
containing the S'-matrix. Then the phase shifts are determined by imposing the continuity of 
x{R) an d of its derivative with respect to R at R = R c . For the NN channel a number of 15 
Gaussians is large enough to obtain convergence and the matching radius is R c = 4.5 fm. The 
size parameter of the Gaussians is fixed at b = 0.44 fm, as discussed above. An important 
analytic step in solving the Eq. (11) is the calculation of two-body matrix elements contained 
in the Hamiltonian kernel 7i(R, R') and the normalization kernel M{R, R'). All relevant two- 
body matrix elements obtained by integration in the color, spin and flavor spaces are given 
in Table 1. Symbolically an upper index / is introduced in order to distinguish between the 
color and the flavor \{ SU (3) matrices. 

Results 

In Figs. 1 and 2 we plot the calculated 3 Si and 1 S'o as a function of Ei a b. For a 



comparison, we add the corresponding results of [^2| obtained with the spin-spin part of 
the one-gluon exchange (OGE) model. The phase shifts of both models reveal the presence 
of a short-range repulsion in the NN interaction. In the GBE model, without the long- 



range Yukawa part the repulsion is stronger and corresponds to a hard core radius 
r^ BE = 0.81 fm (versus r$ GE = 0.35 fm) for and rfi BE = 0.68 fm (versus r$ GE = 0.30 
fm) for 3 S\. This outcome is consistent with the findings of Ref. where a simplified 
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SU(2) version of the GBE model has been used. Thus the GBE model can explain the 
short-range repulsion, as due to a flavour-spin quark-quark interaction and to the quark 
interchange between clusters. However, to describe the scattering data and the deuteron 
properties, intermediate- and long-range attraction potentials are necesary. In calculations 
based on OGE models they were phenomenologically simulated at the NN level by central 
and tensor potentials respectively f23| (see also Ref. ||). However, for a consistent picture it 
is desirable to search for the origin of the attraction at the quark level. Here we incorporate 
a a- meson exchange interaction at the quark level and study its effect on the 1 Sq phase 
shift. Note that the 1 Sq phase shift is not influenced by a tensor potential. For the cr-meson 
exchange interaction we choose the following form 

K = (- ) , (15) 

with parameters discussed below. The introduction of such an interaction is consistent with 



the spirit of the GBE model. It simulates the effect of two correlated pions [|14]]. The good 
quality of the baryon spectrum is not distroyed by the addition of this interaction which 
essentially leads to an overall shift of the spectrum |24j|. The sensitivity of the 1 Sq phase 



2 

shift with respect to the coupling constant -f 3 -, the mass /v and the cut-off parameter A CT 
can be seen from Figs. 3-5. As expected, the attraction in the NN potential increases with 

g 2 q 2 

■j 2 and hence the value of E iab where the phase shift changes sign also increases with as 



shown in Fig. 3. Note that the potential V a of (p~5|) remains attractive as long as /i CT < A CT . 
However fi a cannot be too close to A CT . As suggested by Fig. 4, the attractive pocket in 
the NN potential becomes too small for fi a > 0.65 GeV, making the repulsion dominant 



2 



and leading to negative phase shifts at all energies, when = 1.24 and A a = 0.830 GeV 



A large difference between \i a and A CT is not good either. From Fig. 5 one can see that 

when A CT > 0.95 GeV and fi a = 0.6 GeV an undesired bound state in the 1 S phase shift is 

g 2 

accomodated at -f 2 = 1.24. This is due to the fact that the contribution of the second term 

4-7T 

in the right hand side of ( JI5| ) becomes negligible and V a brings too much attraction in the 
NN potential. 



In this way we found an optimal set of values 



2 „2 

1.24, fM a = 0.600 GeV , A CT = 0.830 GeV . (16) 



9oq 9irq 
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As one can see from Fig. ||, with these values the theoretical curve gets quite close to 
the experimental points without altering the good short-range behaviour, and in particular 
the change of sign of the phase shift at ~ 260 MeV. Thus the addition of a a-meson 
exchange interaction alone leads to a good description of the phase shifts in a large energy 
interval. One can argue that the still existing discrepancy at low energies could possibly be 
removed by the coupling of the 5 Dq N-A channel, suggested by Ref. [^5| in the frame of a 



hybrid model, containing both gluon and meson exchange at the quark level. To achieve this 
coupling, as well as to describe the 3 Si phase shift, the introduction of a tensor interaction 
is necessary. 

Conclusion 

This study shows that the GBE model gives good promise in describing the NN scattering 
properties, while offering a good description of baryon spectra. In the following step it will be 
interesting to include a tensor interaction at the quark level, which does not alter the quality 
of the spectrum [26| and which provides the long-range attraction in the NN potential. 
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FIGURES 
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FIG. 1. The 3 S\ NN scattering phase shift as a function of £7^. The solid line shows the 
result obtained in the GBE model, the dashed line in the OGE model of Ref. [19]. 
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FIG. 2. Same as Fig. 1 but for the ^0 NN partial wave. 
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FIG. 3. The 1 Sq phase shift for various values of g = g£ g /47r at fixed [i a = 0.600 GeV and 
0.830 GeV. 
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FIG. 4. The 1 Sq phase shift for various values of m = n a at fixed g^ q /^TT = 1.24 and A CT = 
0.830 GeV. 
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FIG. 5. The 1 Sq phase shift for various values of I = A a at fixed <j£ g /47r = 1.24 and /i CT 
0.600 GeV. 
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FIG. 6. The 1 5o NN scattering phase shift obtained in the GBE model as a function of E\ a \,. 
The solid line is the same as that of Fig. 1 and the dashed line includes the effect of the u-meson 
exchange potential (12) between quarks with \x a = 0.600 GeV and A CT = 0.830 GeV. Experimental 
data are from Ref. [25]. 
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TABLES 



TABLE I. Matrix elements (NN\0\NN) of different operators O for (S, T) = (1, 0) and (0, 1). 



o 


(5,T) = (1,0) 


(S,T) = (0,1) 


1 


243 


243 


p/cc 


-3 


-3 


AJ • A2 


-648 


-648 


A3 • Ag 








\c \c p/o-c 

Ai • An -Too 
1 ^ OO 


8 


8 


\c \c p/°" c 

O D OD 


-16 


-16 


\c \c p/crc 
A l ' A 3 36 


8 


8 


\c \c p/o"C 
^1 ' ^6 -'36 


8 


8 


\c \c p/cc 

AI • A4 -Tq fi 

^ *± ou 


-4 


-4 


0"! • 02 


-243 


-243 


3 • 06 


27 


-81 


01 • 02 


51 


51 


3 • 6 PL 

' J ou 


57 


93 


01 • 3 P£ 6 


-21 


-21 


01 • 06 -P36 
x ' ou 


-21 


-21 


r>f°" c 

01 • 04 r, fi 
A ^ ou 
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01 - 02 T\ • T2 


1215 


1215 


0"3 ' 06 r 3 ' r 6 


-225 


-225 


nf CTC 
01 • 02 Ti • T2 -Tog 


-111 


-111 


0"3 • 06 T3 • T6 ^36 


177 


177 


01 • 03 n • r 3 FjJq 


33 


33 


01 • 6 Ti ■ r 6 P 3 J 6 


33 


33 


01 • 4 Ti • T 4 


9 


9 


01 • 02 A{ • A^ 


1134 


1134 


3 • 6 A{ • A^ 


-216 


-252 


01 • 02 a{ • a£ p 3 / 6 ctc 


-94 


-94 


03 • 6 A{ • A^ p(£ c 
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208 


\f \ / r>f ac 
01 • 03 A x • A 3 F 36 


26 


26 


01 • 6 A{ • A^ 


26 


26 


01 • 04 a{ • a{ p(q c 


11 


9 


factor 


1 

243 


1 
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